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The SCIP2 Program Documentation is divided into two sections: 
the User’s Guide and the Program Description. Section I details the 
operation of the SCIP2 program with a brief general description and 
purpose of the program. Section II describes the program structure 
in detail. A summary of the equations coded in the program is found 
in the Appendix. An accompanying document, NASA CR-112121, contains 
the technical details of analytical techniques and a preliminary study 
which makes use of these techniques. 
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INTRODUCTION 

SCIP2 is a digital computer program which can be used to investigate the 
effects of instrumentation errors on the accuracy of aircraft stability and 
control derivatives identified from flight test data. The program is based 
on the assumptions that the aircraft differential equations of motion are 
linear and consist of small perturbations about a quasisteady flight condi- 
tion. It is also assumed that a Newton-Raphson optimization technique is used 
for identifying the estimates of the parameters. The full theoretical develop- 
ment of the program equations appears in a companion document entitled "Analysis 
of Instrumentation Error Effects on the Identification Accuracy of Aircraft 
Parameters". A summary of the equations which are coded in the program are 
contained in the Appendix of this User's Guide. 

The uses which can be made of SCIP2 include: 

1. The determination of the effect of instrumentation errors on the 
statistical accuracy of the stability and control derivatives and 
other parameters identified from flight test data can be made. 

This includes the mean error and standard deviation of each of the 
parameters identified. The contribution of each error source on 
each parameter is determined. 

2. The effects of such variables as aircraft type and flight condition, 
control input sequence, and data sampling rate on the accuracy of 
the identified parameters can be determined. 

3. Tradeoff studies can be made between instrument quality and 
identification accuracy. 

4. Different combinations of instruments can be studied for use in 
collecting the flight data. 

5. Tradeoff studies between fewer instruments with greater quality and 
more instruments can be made. 
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The necessary instrument accuracy required in a flight test program 
to allow identifying aircraft parameters to a desired level of cer- 
tainty can be specified. 


PROGRAM EXECUTION INFORMATION 


2 . 1 Machine Requirements 

This program was developed on the UNIVAC 1108 machine and then converted 
to the CDC 6600. Thus, it may be executed on either machine with very little 
change. 

2.2 Storage Requirements 

It requires (71,500)g words to load and (56,600)g words to execute this 
program on the CDC 6600. This requirement results from compiling the program 
under FORTRAN 2.3 RUN COMPILER. The core requirement may be reduced by re- 
moving double precision arrays or by assigning a single buffer to two scratch 
files. 

It is also possible to reduce core requirement by using the overlay 
feature of the SCOPE operating system and also by using blank commons in place 
of labeled commons. 

2.3 Execution Time 

The execution time for the program varies with the various options of 
the program. The shortest case will be the ESTIM=1 Short Period Ensemble 
Analysis without lag, and the largest will be the ESTIM=3 Lateral Simulated 
Data Analysis with lag. 

It took 6.55 seconds of CPU time to execute a five second trajectory Short 
Period ESTIM=4 Ensemble Analysis without lag, whereas 14.28 seconds to 
execute the same case under M0DE=1 of the Simulated Data Analysis. 

The execution time increases if lag is present since it requires to 
integrate equations with a smaller integration step. 


2. 4 Multiple Case 


The program allows execution of more than one case per job. This is 
accomplished by having sets of input data one after another. Each case 
starts with POD card (see Section 3.0). 
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INPUT ELEMENT DESCRIPTIONS 

There are three methods for providing input data to the SCIP2 program: 

1. By means of the BLOCK DATA subprogram, 

2. By data card input at execution time, and 

3. By a combination of 1. and 2. above. 

3 . 1 Data Card Input 

Data card input is optional except for the program options data (POD) 
card. All data not set in the BLOCK DATA subprogram or data to be changed 
may be flagged for input on the POD card. The required POD card has a fixed 
field format as follows: (A6,4X,30L1) . 

The POD card selects the error analysis type: ensemble (ENSMBL) or sim- 
ulated data (SIMDAT) . Card columns 11 through 40 are divided into 30 logical 
fields of length 1. If a T appears within one of these fields, then the option 
corresponding to that field is exercised. If the field is blank or contains 
an F, then that option is ignored. Fields 2 through 22 (card columns 11-32) 
are reserved for input options. Fields 23 through 30 are reserved for other 
program options. 

Table 3.1 shows the correspondence between the input option field and 
the group of data controlled by the option. 

3 . 2 Input Elements 

Table 3.1 lists and describes the input groups (namelist identifier, if 
applicable), variables in each group, the unit, and the type of each variable. 


"Card Field" gives the column number of the input option on the POD card. 
Blank under "Namelist Identifier" indicates fixed field data. Under "Type", 
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the variable types are given as: 


A - hollerith data read in A6 format 
F - floating point number 
I - integer number 

Note that all hollerith data are read in fixed format and all numbers 
are in namelist format. 


3.3 BLOCK DATA SUBROUTINE 


The BLOCK DATA subroutine may be thought of as the SCIP2 data base and 
is useful when most parameters for a study remain fixed for many runs while 
several parameters change for each run of the study. If it is desired to study 
the effects of different instrument groups for a specific aircraft then a BLOCK 
DATA subroutine may be constructed consisting of the nominal aircraft and instru- 
mentation model. Then, for succeeding runs those parameters requiring changes 
may be changed by data card at execution time. 

A nominal BLOCK DATA subroutine is provided with SCIP2. It consists of 
the LABEL COMMON blocks which will contain program inputs. These are defined 
in the Table 3.2. 


NOTE: The following order for parameters and instruments is used by the SCIP2 


program. 
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3.4 Sample Input 
3.4.1 Example No . 1 

The following input is for running the Ensemble Analysis subprogram with 
the full longitudinal equations of motion. Only the parameters of the F and 
G matrices are estimated. Seven instruments are used to measure the output. 
The output from this example is presented in Section 4.2. 


ENSMBL . TT TTTTT TTTTTTT T TT TTTTTT 
FULL 

PAG _ PPG aoav_ avpt _axacc norac paacc_ 

' SRFTRAJ THTet A0=-~.556, 'TLpH A 0 =- . 7, VEL = 178.9, $' 

_ S OP TO NS ESJ IM=1 $ 

$MIS3~ ' ' NMC =1 , D T =9 .05, $' 

SL0N3DR _ __ _ 

MQ=-0. 037694, MW=- 0.05 0 715', ’ Z W = -2 . 07 0 17 4 , MU=-0, 0 00620, 

ZU=-. 334407, XU=-. 042662, XW=0. 070207, M0E=-24.38, ZDE=0 

Tend ’ 

SCONSEQ 

NO? s IOO, 

DELE=2Q*1. 0 ,20*-l. 0 , 100*0.0, 

‘“SEND' 

_$PITATT WTHETA = .5_, 3T HEJA= 0 . 2 , ETH ETA= . 005 $ 

SPf TR AT W Q= . 5 , B Q = 0 .2 , EQ= .00 5“ $" " 

SAOAVAN HAL?HjA=«5_» 3A L_PHA = . 3 , E A L PHA =.005 $ 

‘ $ AVPT US WU = . 5 , " BU-1. , EU= . 00 5 $ 

SAX I ACC WNXf.l, 3 NX = . QQ1, ENX = . 0 0 5,_G N X=. 5$ 

$NO RM AC~ WNZ= . 0 25 , "BN Z = .001 , E NZ= . 005, GN Z = . 5 $ 

SPANACC WQDDT=.5 ? _ 3Q9 0^= • 2 L Q00 I = • 0 0 5 -_?- 
” $CG UN C R E XC j = • 5", EZCG=.5~$ 

3CG A 3 L E A X=2 • , EAZ=2. 3 
$ A 3 VANE' EV X= 5 T" $"" 


3.4.2 Example No. 2 

The following input is for running the Ensemble Analysis subprogram with 
the lateral equations of motion. The parameters, the initial conditions, and 
all the output instrument biases except b^ are estimated. Seven instruments 
are used to measure the output. 


ENSMBL TT T TTTT TTTTTTT T TT TTTT TT ... 

LAT 

!AOSV__ RRG YRG R A G LA T A C RAA CC YA ACC 

"SRFTRAJ THET A 0=- • 556 , ALPHA 0=-.7, VEL=178.9, $ 

SOPTQNS ESriM =4 , $ 

S Ml S3* D T = 0.1, $ 

_3LATDR . .. 

YBETAa-,163, LOET A=-23 .26 , N3ETA=5.504, LP=-il.53, NP=-i.363, 

LR=2 .69 2, NR=1 . 214, Y O ELA= Q. , L0 ELA=53.79, ND_ELA_=.2103_, N0ELR=-6. 17 2 
' XZ 1= 0. , ‘ 

SEND ; • _ 

“ SCONS EQ ‘ * ’ * * 

MOP_=50 l 

DELR=2 0*0. 5, 5*- 1.0, 5*1.0,100*0. 0, 

D EL A.=5 5 0 .J j 5_*-3_._5 , 20*-0 .0 5, 10 0*0. 0 , 

'SEND ’ 

_ SAOSy AN_ W9ET A = . 0 09, 8CET A = 0. , E3£TA=.0Q5, S 

SROLRAT WP =.00*5, BP =0^ , EP =.005 , GP =.5, $ 

SYAWRAT _WR . =.005 , 9 R =0. , ER =. 005, GR =.5, $ 

~ SROLATT WPHI =.0 05, BPHI =.2, EPHI =. 035 , $ 

3 LATA CC WNY =._0_0 1 , . 3NY =0 . , ENY_. =. 0 05 , $ - 

SRANACC WPDO T= .0 05, BPDOT=0. , EPOOT=. 005 , GP0OT=.5, $ 

SYAJNA CC __W_RD OT=^0 05, _8 RO 0 T =0., ERDOT=.OQ5 j G R_DOT = . 5*. 

'SCGUNCR E XC G = • 5 , "EZCG~ =• 5 , $ 

_ SCS AP L EAX f 2 • 0,_ E_AZ = 2.J1,__ $ 

$ A3 VANE EVX =5.0, $ 


3.4.3 Example No . 3 

The following input is for running the Ensemble Analysis subprogram with 
the longitudinal equations in the short period mode. Only the parameters are 
estimated. Four instruments are used to measure the output. The angle-of- 
attack vane is omitted. 

ENSMOl TT TT TTT T T T T T T TTTTTT 

SPL " ' ““ 

PAG ,._PRG NORA C_PA ACC 

SRFT.RAJ THETA 0=- .556 , ALPHA0 = -.7, V £1=178. 9, $ 

SOPTDNS ESTIM=1, $ 

SMISS NMC=1 , D T = 0 . 0 5 , $ 

SLON3DR 

HQ=-0. 037594, MW=- 0.050 715 ‘Tw=-2 ^ 07 0 1 747 ~ 

M0E=-24. 330927, ZDE=0.0, 

SEND : 

SCONS EQ 

D EL E= 20*1.0,20^-1.0, 100*0.0, ' “ 

SEND 

SPITATT WTH ETA = 0 • 25 , 3THETA = 0'.l, E tWT A = 6 7605*J $ 

_S _P I T J _ W Q=0.?,_BQ=g. 1, E Q = 0.0 05, _$ 

$ NORM AC WNZ =0.025, BNZ = 0.0*i, E NZ= 6 .005, G NZ =0.57 $~ 

_IPANACC_ WQOOT = 0. 3, 3Q.QOT=0.1, EOnOT = 0,Q05, $ 

SCGOnCR EX3 G = 0 . 5 , '$. 7 

SCGA 3 L E AX= 1 . Q , $ 


3.4.4 Example No. 4 


The following input is for running the Simulated Data Analysis subprogram 
with the longitudinal equations in the short period mode. This is a special 
case with no instrument errors. Two runs are made with parameters perturbed # 

5% in the second run. The output is the result of applying the modified 
Newton-Raphson algorithm to estimate the difference in the parameters. 


S I HD AT ft TTTTT TTT ft T TT TTTT "* • 

SPL _ 

PAG PRG A 0 A V NO RAC PA ACC ' 

SRFTR A J THETA0=-.55£> , ALPHA0=-.7, VEL=178.9, $ 

$ OPT 5 NS’ ESfl H=i, J 

SMI SC NMC = 1 i OT=0.35, MOOE=2, $ 

$ LONG DR “ * ' “ • 

MQ=-0. 037594, MW=- 0 . 350 715 , Z W=-2 . 07 0 1 74 , 

MDE =-24. 38 0 927", ZQE^OVO, “ 

_$END __ 

SCONS FQ " “ ’ 

DELE = 20*1. 0 ,20*-l. 0 , 100*0.0, 

'"SEND "" ‘ ' ~ • 

SPITATT WH E T A=0 . 25 » $ 

$ P I TR AT ~W Q= 0 • 2 f ' $ ~ 

SA0A7AN _WALPH A=0_.5_, _$ 

S NO PM AC HNZ= 0.02 5 , $ 

SPANACC W GO 0 T = 0 . 3 , $ 

SI MOAT ft f T TTTTTT * • 

SPL _ _ __ 

PAG PRG ' A 0 A V ' NOR AC PA ACC * ' 

SLONG DR 

MQ= - 0 . 0 395 , MW=-0. 0533, ZM=-2.18, 

MD E = - 25 . 3 , ZQE = 0 «_0 , 

SEND ~ ” • 


3.4.5 Examp le No . 5 


The following input is for running the Simulated Data Analysis subprogram 
with the short period (longitudinal) equations. This is an example of Mode 1, 
and each of the error sources is applied individually to the output measure- 
ments. The contribution due to noise is set equal to 
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S I MOAT TT TTTTT TTT TT T T TTTTTT 

SPL 

PAG PRO AOAV . : NORAC PAACC 

"fRFtRAJ THET A0=- • 556 , ALPHA 0=-. 7, VEL=178.9, $ 

$ 0 P T 0 NS__ tSTl M~l> $ : 

SMISC NMC=1, DT=0.05, MODE=i, $ 

S L 0 NO OR __ 

MQ=-0. 037594, MW=- 0.350 715 , Z W=-2 . Q70 174, 

M 0 £ = - 2 4.33092 7, Z DE= 0. 0 , 

SEND 

JCON3EQ . 

DEL E=2 6* 1. Q»20*-1.0 , 10 0*0.0, 

SEND 

-jpffATT WT-j E T A=0 • 25 , BTHET A=Q .1 , ETHET A = 0. 005, $ 

SPITRAT _WQ=0_.?,_ Bq=_0._l,_ EQ = 0.0 0_5_, $ 

SAOAVAN WALPHfl = 0.5, DALPHA= 0 . 1 , EALPHA=0. 0 05, $ 

' SNORHAC WNZ =0.025, 8NZ = 0_.0_1 , __EN7=0 . 005 , . GNZ=0.5, S 

SPANACC WQDOT=0.3, BQDOt=0. 1, EQDOT=Q. 005, $ 

3CGUNC R E XCG=0.5 ♦ ? 

'"SCGAPL t ax =1 . o r~ ' 



The following input is for running the Simulated Data Analysis subprogram 
with the short period equations. This is an example of Mode 2, which is the 
Monte Carlo option. Twenty Monte Carlo runs are made. Mean errors are constant 
in each run. Random errors except noise are randomly set at the beginning of 
each run and held constant during the run. Noise errors are randomly set at 
each sample point. 


SIMOAf TT TTTTT TTT TT T TT TTTTTT 

■"SPL'' “ 

PAG _PRG AOAV NORAC PAACC 

SRFTRAJ THET A0=-".556 , ALPHA 0 = -. 7, VEL=i78.9, $ ~ 

SOPT3NS ESTJM=1, _S 

• SMISC NMC=20, M~ODE=2, DT=0.05i $ * 

SL ON3 or 

MQ=-0 . 037594, MW=- 0 . 0 50*715 , Yw=-2~. 070174, 

MDE=-24.330 927, ZDE=0.0, 

send' * " " “ 

_SC0N3EQ 

• DELE=20*1. 0,20*-l/0, 100*0.0, ‘ ' 

SEND 

SPITATf W THE T A = 0 • 2 5» - 9THET A~=0 • 1 , E THET A = 0 .0 0 5 , $ 

__$P IT R A T WQ= 0. 2, • BQ=0.1, EQ = 0. 0 0 5 , $ 

SAOAVAN WALPH A=0 • 5 , DALPHA=d.l, E ALPHA = 0.0 6 5, $ 

$NORM_AC WN7=0. Q25_, 9NZ = 0_._0_1 , EN7=0. 005, _GNZ=0.5, S 

• SPANACC WQDbT = 6.3, BQQOT=0.1, EQDOT=0.005, $" 

SCGUNCR EXCG=0.5, $ 

$CGA 3 'l ' EAX=1.0, $ 

SA P VAN E EVX= 1 0.0, S 
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3.4.7 



The following input is for running the Simulated Data Analysis subprogram 
with the lateral equations. This is an example of Mode 3, which is just one 
pass through the program with a fixed set of instrument errors. 


si MO AT TT ft iff” T T TT TTT T TT TTTTTT 
LAT 

AOSV ~RRG Y RG RAG~' LAT AC RAACC YAACC ' ’ ' 

SRFTRAJ THET AC =- , 5 56 , AL PHA0--.7, VEL=178.9, S _ 

SOP TONS EST IM=i, $ 

jot=o.i, mode= 3| s ______ 

SLATDR “ * ' 

YBET A = -. 163 , LJ3ET A=j-23L.25 , N3 ETA = 5_. 50 4 , LP=-ii.53, NP=-1.363, 

LR=2.692, MR=1.214, YDELA=0. , LDELA=53. 79, NOELA=.2103, NDELR=-6. 172 , 

XZI=_0. t _ __ __ 

SEND 

_$CONSEQ_ __ 


NO P= 50 

* 




D EL R= 

20*0.5,5*-1.0 

, 5*1.0,10 0*0. 0, 



DEL A = 
SEND 

5*0.5 ,5*-1.0, 

2 0*-0 • 05, 10 0*0.0, 



SA0S7 AN 

WBETA=.0Q5, 

BEET A = . 3 , E3ET A= . 005 , 

$ 


SROLRAT 

WP =*005, 

8P =.2, EP =.005, 

GP =.5, 

S 

$YA WRAT 

WR =.005, 

BP =.2, ER =.005, 

GR =.5, 

S 

SROL4 TT 

WPHI =. 0 05, 

BP9I =.2, EPHI =.005, 

$ 


SLA T A CC 

WNY =.001, 

9NY =.001, ENY =.00 

5, $ 


sranacc 

WP3 0 T = . 0 05 , 

BPD 0 T = . 2 , EPDOT=. 005 , 

GPDoT= .5, 

S 

SYANACC 

WR3QT = .0 05, 

OROOT=.2, ERDOT=.005, 

GRD0T= .5, 

$ 

SCGUNCR 

EXCG =.5 , 

E7CG =.5 , $ 



SCG A 3 L 

E AX =2.0, 

E AZ = 2. 0, S 



SA3VA NE 

EVX =5.0, 

S 




3.5 Discussion of Input Options 

Several of the input options require elaboration. Each will be discussed 
below, separately. 


12 


3.5.1 Inputs (23) = T 

This is the special input option. By setting INPUTS (23) = T (field 33 
on the POD card), the user may iterate through the Simulated Data Analysis 
portion of the SCIP2 program. In the iterate mode, the program makes a single 
pass through the ensemble analysis (retaining the required partial derivatives 
and covariance for the simulated data pass). From this point on, the SCIP2 
program processes the new data only in the simulated data analysis and only 
for a single iteration, always using the necessary results from the ensemble 
analysis generated from the original input data. This allows the user to 
vary any of the input parameters from their nominal values and study the effects 
via the simulated data analysis. (Note that, if any of the aircraft para- 
meters are changed on subsequent iterations all must be re-input because of 
units conversions within the program. All other inputs carry no added restric- 
tions in this iterate mode.) 

In particular, the aircraft parameters may be perturbed from their 
nominal values in this wav and the effects studied (see Section 3.4.4, 

Example No. 4). 

3.5.2 PERT and NDXP 

If INPUTS (23) = T, (that is, the program is being operated in the 
iterate mode), then the aircraft parameters in the Lateral equations may be 
perturbed from the nominal input by the percentage indicated by PERT. If 
NDXP = 0, then all lateral parameters are perturbed one at a time (the others 
fixed at nominal values). If NDXP = 1 <_ n <_ 13, then the n 1 "* 1 lateral para- 
meter only is perturbed. 

3.5.3 NMC = 1, MODE = 1 

If NMC = 1 and MODE = 1 is input to SCIP2, then the covariance matrix 
generated by the ensemble analysis phase is used in the simulated data analysis 
rather than regenerating the covariance matrix using a Monte Carlo process. 

This saves computational time since, theoretically, these covariances are 
identical and since in MODE 1 it is desired to examine the effects of the 
error sources other than white noise. 
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Table 3.2 


LABEL COMMON 
BLOCK NAME 

VARIABLE OR 
ARRAY NAME 

DIMENSION (IF 
APPLICABLE) 

DESCRIPTION 

BIASES 

BLONG 

7 

longitudinal instrument biases 


BLAT 

7 

lateral instrument biases 

CGNCRT 

ECGVEC 

3 

C. G. uncertainties (e e e . 

cgx, cgy , cgz) 

CNFLCT 

SPLCNF 

3 

short period initial conditions conflict 
vector. SPLCNF (n)=T implies that the 
nth initial condition correlates with 
a bias error and is to be treated as 
an error source. 


FULCNF 

4 

full longitudinal conflict vector, 

FULCNF (n)=T implies that the nth initial 
condition correlates with a bias and 
is to be treated as an error source. 


L AT CNF 

4 

same as FULCNF but for lateral analysis 


SPLBCN 

5 

short period bias conflict vector, 

SPLBCN (n)=T implies that the n fc h bias 
correlates with an initial condition 
and is to be treated as an error source. 


FULBCN 

7 

Same as SPLBCN but for the full longi- 
tudinal analysis 


LATBCN 

7 

same as SPLBCN but for the lateral 
analysis 

CONTRL 

NP 

1 

number of data points to the end of 
trajectory at intervals of DT 


DELE 

300 

NP elevator deflection at intervals of 
DT 


DELA 

300 

NP aileron deflections at intervals 
of DT 


DELR 

300 

NP rudder deflections at intervals 
of DT 

CPERRV 

WDE 

1 

Standard deviation for control measure- 
ment noise for 6 

e 


WDA 

1 

Similar to WDE but for 6 

a 


[ 
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Table 3.2 (Continued) 


LABEL COMMON 

VARIABLE OR 

DIMENSION (IF 



BLOCK NAME 

ARRAY NAME 

APPLICABLE) 

DESCRIPTION 

• 


WDR 

1 

Similar to WDE but for 6 

r 



BDE 

1 

Bias error for control measurement 
for 6 





e 

• 


BDA 

1 

Similar to BDE but for 6 

a 



BDR 

1 

Similar to BDE but for 6 

r 

• 


EDE 

1 

Scaling error for control measurement 





for <5 

e 



EDA 

1 

Similar to EDE but for 6 

a 

• 


EDR 

1 

Similar to EDE but for 6 

r 



FDE 

1 

Control measurement lag for 



FDA 

1 

Similar to FDE but for 6 

a 

• 


FDR 

1 

Similar to FDE but for 6 

r 


DDDPBK 

ROWDLO 

9 

row numbers for which a non-zero 





element of 9D/9p exists (longitudinal 
analysis) 

• 


ROWDLA 

13 

row numbers for which a non-zero 





element of 9D/9p exists (lateral 
analysis) 

• 


COLDLA 

13 

same as ROWDLA but for the column 
number 


DFDPBK 

ROWFLO 

9 

row numbers for which non-zero element 
of 9F/9p exist (longitudinal analysis) 

<9 


COLFLO 

9 

same as ROWFLO but for column numbers 



ROWFLA 

13 

row numbers for which non-zero element 
9F/9p exist (lateral analysis) 



COLFLA 

13 

same as ROWFLA but for column numbers 

• 
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Table 3.2 (Continued) 



LABEL COMMON 

VARIABLE OR 

DIMENSION (IF 



BLOCK NAME 

ARRAY NAME 

APPLICABLE) 

DESCRIPTION 


DGDPBK 

ROWGLO 

9 

row numbers for which non-zero element 
of 3G/3p exist (longitudinal analysis) 



ROWGLA 

13 

row numbers for which non-zero element 





of 3G/3p exist (lateral analysis) 



COLGLA 

13 

same as ROWGLA but for column numbers 


DHDPBK 

ROWHLO 

9 

row numbers for which non-zero element 





of 3H/3p exist (longitudinal analysis) 



COLHLO 

9 

same as ROWHLO but for column numbers 



ROWHLA 

13 

row numbers for which non-zero element 
of 3H/3p exist (lateral analysis) 



COLHLA 

13 

same as ROWHLA but for column numbers 


DMATRX 

DLONG 


D-matrix for the longitudinal 
analysis 



DLAT 


D-matrix for the lateral 





analysis 


FMATRX 

FFLONG 

4x4 

F-matrix for the longitudinal analysis 



FFLAT 

4x4 

F-matrix for the lateral analysis 


GMATRX 

GLONG 

4 

G-matrix for the longitudinal analysis 



GLAT 

4x2 

G-matrix for the lateral analysis 



XXI 

1 

Moment of inertia I 

XX 



XZI 

1 

Product of inertia I 

xz 



ZZI 

1 

Moment of inertia I 

zz 


HMATRX 

HLONG 

7x4 

H-matrix for the longitudinal analysis 



HLAT 

7x4 

H-matrix for the lateral analysis 


LAGERR 

FLONG 


diagonal elements of the lag 

error matrix F for the longitudinal 

analysis 

• 


FLAT 


diagonal elements of the lag error 

matrix F for the lateral analysis 
m 
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Table 3.2 (Continued) 


LABEL COMMON 
BLOCK NAME 

LOCCG 

MISALN 


MISCEL 


MISCL1 


VARIABLE OR DIMENSION (IF 

ARRAY NAME APPLICABLE) DESCRIPTION 


CGLOC 

EVX 

GMLONG 


GMLAT 


DT 

NMC 

GRAV 

UNITS 

RADDEG 

NDXP 


accelerometer location error vector 

(e , e , e ) 
cgx * cgy cgz 


vane location error 

misalignment error vector ( Y 0 , , 

Y , y > Y » Y > Y* ) for the 
a * u nx 'nz’ 'q 

longitudinal analysis 

misalignment error vector (y„, y , 

p P 

Y , Y. » Y > Y * j Y * ) for the lateral 
r <{> ny p r 

analysis 

integration step size 

number of Monte Carlo runs to be 
executed (if ENORSD=SIMDAT) 

gravitational constant 

units of gravitational constant 

radians per degree 

index for perturbing lateral 
parameters in iterate mode 


PERT 

TIME 

DTS 

TAU 

NUM1 


percentage perturbation of lateral 
parameters in iterate mode 

time of the sample point 

time between the two successive sample 
points 

time to the first sample point for 0 

lagged controls 

number of elements in the control 
variable 


NUM2 


not used 
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Table 3.2 (Continued) 


» 


LABEL COMMON VARIABLE OR DIMENSION (IF 

BLOCK NAME ARRAY NAME APPLICABLE) DESCRIPTION 

MCDM2 variable MODE - 2 

II running counter of the Monte Carlo 

runs 


SWICH 


IU1 


flag, if non-zero, to indicate that 
the initial condition eror is being 
considered 


scratch unit to save 



y» 



ID 2 


scratch unit to save 6p for each of 
MODE = 2 Monte Carlo runs 


NAMVAR NAMVAR 

NAMIC 

NAMVAN 

NAMCG 

NOISE WLONG 

WLAT 

OPTNS ENORSD 

EQNS 


7x2 


4x2 


4 


3 


7 

7 


an array containing the hollerith names 
of the measured output states for longi- 
tudinal and lateral analysis 

an array containing the hollerith names 
of the aircraft states for longitudinal 

hollerith names 
e 

vx 

hollerith names 
fe e 
' cgx * cgy' 

noise vector 

noise vector 

for ensemble 
analysis, SIMDAT for simulated data 
analysis 

equations option - FULL for full 
longitudinal analysis; SPL for short 
period longitudinal analysis; LAT for 
lateral analysis 


and lateral analysis 

an array containing the 

for (e e , e ) and 
ax’ ay’ az 


an array containing the 

for C-G location error, 

e ) 
cgz 


standard deviations for 
(longitudinal) 
standard deviations for 
(lateral) 

program option - ENSMBL 
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Table 3.2 (Continued) 


LABEL COMMON VARIABLE OR DIMENSION (IF 

BLOCK NAME ARRAY NAME APPLICABLE) DESCRIPTION 


ESTIM 


parameter estimation option: 

1. estimate parameters only 

2. estimate initial conditions also 

3. estimate biases as well as 1 & 2* 

4. estimate i.c.'s as well as 1 and 

biases* 


MODE - if ENORSD = SIMDAT, then MODE represents 

the Monte Carlo options: 

1. Monte Carlo on noise 

2. Monte Carlo on noise and random 
errors 

3. one pass with all error simultaneously 


INSLON 


INSLAT 


7 


7 


desired instruments (longitudinal) 


INSLON 

(1) 

= PAG 

INSLON 

(2) 

= PRG 

INSLON 

(3) 

= AOAV 

INSLON 

(4) 

= AVPT 

INSLON 

(5) 

= AXACC 

INSLON 

(6) 

= NORAC 

INSLON 

(7) 

= PAACC 

desired instruments (lateral) 

INSLAT 

(1) 

= AOSV 

INSLAT 

(2) 

= RRG 

INSLAT 

(3) 

= YRG 

INSLAT 

(4) 

= RAG 

IN SLAT 

(5) 

= LATACC 

INSLAT 

(6) 

= RAACC 

INSLAT 

(7) 

= YAACC 


RANDOM IRN - random number generator starting 

number (UNIVAC-1108 only) 

* Estimate option 3 estimates those i.c.'s which do not correlate with any bias. 
Estimate option 4 estimates those biases which do not correlate with any i.c.'s. 
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Table 3.2 (Continued) 


LABEL COMMON 
BLOCK NAME 

VARIABLE OR 
ARRAY NAME 

DIMENSION (IF 
APPLICABLE) 

DESCRIPTION 


REFTRJ 

VELOC 

3 

reference trajectory 

vector (a , 0 , V ) 
o o 

parameter 

SCALING 

EPLONG 

7 

^ e 0* e q* V e u’ e nx* 

e nz’ e q^ 


EPLAT 

7 

(eg, e p , e r , e^, e ny . 
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OUTPUT DESCRIPTION 


4.1 Printed Output Control 

The output from this program is controlled by elements 25 through 30 of 
the INPUTS vector. If the element is "TRUE", then the output associated with 
it is produced, otherwise not. The output quantities of both the Longitudinal 
or Lateral analysis are the same except for the parameter headings. 

4.1.1 Measurement Vector 

If INPUTS (25) is "TRUE", then the measurement vector y is printed at 
every sample point (every DTS-seconds) . This applies to both subprograms 
(e.g., ENSEMBLE ANALYSIS or SIMULATED DATA ANALYSIS). 

4.1.2 Ensemble Analysis 

Besides the measurement vector y , the following quantities may be 
printed from the ensemble analysis: 


(a) INPUTS (26) controls the printout of the covariance 
matrices. 


a 2 j 


8P 


-1 


and 



5 2 J 

a P" 


^ is the convariance due to noise and 


/ aP \_/ T. / aP \ T is the covariance due to the random 

(te-r 66 ) ^ 


error sources. 


INPUTS (27) controls the printout of the normalized cov- 
ariance matrices described above. The relation between 
the normalized and unnormalized covariance matrix is 
given by following equations: 


diagonal elements C. . = sjc' . 7 


« off-diagonal elements C . . = ij 

1J C:. . C . . 
n jj 


where the primed variables are the unnormalized covariance 
matrix elements and the unprimed variables are normalized. 


(c) INPUTS (28) controls the printout of the partial deriva- 
tive matrix, 3P 
_9e_ 


(d) INPUTS (29) controls the printout of the mean expended 

errors, 3P E(e). 

9e 


4.1.3 Simulated Data Analysis, Mode = 1 


(a) INPUTS (26) controls the printout of the covariance 
matrices , 


3 2 J -1 and 


i i ( p ^ T ) + i(v/) 


where P. is the error vector due to the noise. 

N is the number of Monte Carlo runs. 

— til 

P. is the error vector due to the i of q random 
J 

error sources . 
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If N=l, then 1^ ^ ^P re Pl ace d by 

(h) INPUTS (27) controls the printout of the normalised covar- 
iance matrices described in 4.1.3 (a). 



(c) INPUTS (28) controls the printout of the error vector 

1-1 


3P^ 


(«)' 


esulting from each of the random and mean 


error sources. Each error vector is printed columnwise up 
to 10 error sources across the page. 

(d) INPUTS (29) controls the printout of the mean error vector, 



4.1.4 Simulated Data Analysis, Mode = 2 


(a) INPUTS (26) controls the printout of the covariance matrix, 

i N - - T 

^ Z (P. - P)(P, - P) * 

N i=l 1 

where N is the number of Monte Carlo runs and P^ is the 
error vector resulting from applying all the error sources 
at once, and P is the average error vector described in (c) . 

(b) INPUTS (27) controls printout of the normalized covariance 
matrix described in 4.1.4 (a). 

(c) INPUTS (30) controls the printout of the average error 

N 

vector , 1^ 'ST P . , 

nZ 1 

i=l 

where N is the number of Monte Carlo runs and P^ is the error 
vector resulting from applying all the error sources at once. 


4.1.5 Simulated Data Analysis, Mode = 3 

The error vector is always printed when M0DE=3 of the Simulated 
Data subprogram is executed. 
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4 . 2 SAMPLE OUTPUT 


The following is an example of the output which results when input, such 
as presented in Section 3.4.1, is used. Explanatory comments appear in 
parentheses. 
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THE STATE MEASUREMENT AND CONTROL SEQUENCE WHICH ARE USED TO 
COMPUTE THE INSTRUMENT OUTPUT AND PARAMETER SENSITIVITIES 
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The matrix 
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DETERMINANT OF D2J0P2 o 3.59220404^*08 



V 


PROGRAM FUNCTIONAL ORGANIZATION 

The program consists of two subprograms: the Ensemble Analysis Program 
and the Simulated Data Analysis Program. 


The main program for the Ensemble Analysis is SUBROUTINE ENSMB. It 


evaluates ensemble equations to construct j — j) matrix and then evaluates 
^ 9 (6p) j _ \ 3p / 


The Simulated Data Analysis Program is further subdivided into M0DE1 
and M0DE2 subprograms, SM0DE1 and SM0DE2, respectively. The Simulated Data 

Analysis Program also requires the I — j I matrix, and this task is accomplished 

\ 3 P / 

by executing the ENSMB subprogram first. 


The figures 5-1, 5-2, and 5-3 show the flow of ENSMB, SM0DE1, and 
SM0DE2, respectively. 
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FIGURE 5-1 (Continued) 
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FIGURE 5.2 (Continued) 
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VI 


SUBPROGRAM DESCRIPTIONS 

The subroutines in the program are described in this section. The 
function and argument list for each subroutine are given. 

6.1 BLOCK DATA 

The BLOCK DATA subroutine is used to initialize variables in labeled 
commons at load time (see Table 3.2). This subroutine also includes all the 
labeled commons for the entire program for the storage allocation purpose 
(CDC-6600) . 

6.2 CNSTRT 

This subroutine computes the T-matrix during the Simulated Data Mode=l 
Analysis when random error sources are applied one at a time. 

Calling Sequence - CALL CNSTRT (INDEX) 
where 

INDEX = n means that the n^ random error source is being considered. 
The error sources are numbered as follows: 

1 £ INDEX £ 7 — — >bias errors (b's) 

8 £ INDEX £ 14 => scaling errors (e's) 

15 £ INDEX £ 21 > misalignment errors (y's) 

22 < INDEX < 25 > initial condition error (x 's) 

— — o 

26 £ INDEX £ 28 > center of gravity position errors 

(e ’s) 
eg 

6.3 CONSTT 

This subroutine computes the T-matrix for the Simulated Data Mode=2 
Analysis by applying the errors all at once. 
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where 

XMAT is the M x L matrix of the function 
T is the integration variable, time 
DX is the derivative being computed 
K is not used 

L is the number of columns in XMAT and DX 
M is the number of rows in XMAT and DX 
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6.5 SUBROUTINE DERIVU 


This subroutine evaluates time derivatives of the control variable when 
control output lag is present. 

Calling Sequence - CALL DERIVU (UL, ULDOT) 


where 


UL is the control vector, u^. 

ULDOT is the time derivative of the control vector, (u^) 


Equation 


“l - - F c “L + F C U I 


6.6 SUBROUTINE DERVU1 

This subroutine evaluates time derivatives of the control vector as well 
as those of the partial derivatives of the control vector with respect to 
error sources. 

Calling Sequence - CALL DERVUl(UL, ULDOT, PDUDE , PDUDED , NX) 


where 

UL is the control vector, u^. 

ULDOT is the time derivative of the control vector, ~ (u T ) 

at L 


PDUDE is the partial derivatives of the control vector with 
respect to error sources, 

PDUDED is the time derivative of PDUDE, 


Co Ul III 

••&) ■ 


j-W 

* dt \9e ) 


NX is the number of components in the control vector, u . 

Li 


Equation: 


\ ' - F C U L + F c U I 


’V 


dt \3e 


3 u t 3F_ 

*cir-+ir ( ' u l + u) 
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6.7 SUBROUTINE DERIVX 


This subroutine computes the time derivative of the aircraft state 
vector, x. 

Calling Sequence - CALL DERIVX (X, XDOT) 

Where 

X is the aircraft state vector 

XDOT is the time derivative of the aircraft state vector 

Equation: 

x = Fx + Gu 

6.8 SUBROUTINE DERIVY 

This subroutine evaluates the time derivative of the lagged output 
measurement state vector, y^. 

Calling Sequence - CALL DERIVY(YL,YLDOT,YI,X,U) 
where 

YL is the lagged output measurement state vector, y . 

Li 

YLDOT is the time derivative of y L> 

YI is the indicated output vector, y^. 

X is the aircraft state vector, x. 

U is the input control vector, u. 

Equation: 

" ' F M y L + F M y I 

where 

y x = T(Hx + Du) 
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6.9 DERIV1 


A 

This subroutine computes the time derivative of the state vector, x: 

4r (x) = Fx + Gu . 
dt m 


Calling Sequence - CALL DERIV1(F,X,G,U,L,M,N,DXDT) 


where 

F is the matrix F 

X is the state vector x 

G is the matrix G 

U is the control vector u 

m 

L, M, and N are used to specify sub-matrix of F and G which 
are used; i.e., L x M sub-matrix of F and L x N 
sub-matrix of G. 

A J A 

DXDT is the time derivative of state vector x, -r (x) . 

dt 

6.10 DERIV2 

This subroutine computes the time derivatives of the partial derivatives 
of the state vector with respect to the parameters 



Calling Sequence - CALL DERIV2(FMAT,PDXDP,X,U,ROWl,COLl,ROW2,COL2,M,N, 

L,DXDPDT,NEQ) 


where 


FMAT is the matrix F 

c^X 

PDXDP is the matrix of the partial derivatives, -r— . 

dp 

A 

X is the state vector, x 
U is the control vector, u 

m 


fl 


( 3F 

matrix, i.e., ROWl(i) = m and COLl(i) = n 


ROW 2 and COL2 are similarly used to specify the non-zero element 
of matrix. 

M, N, and L are used to specify the sub-matrices of F and G 
which are used; i.e., M x N sub-matrix of F and 
M by L sub-matrix of G, respectively. 

d / 3x\ 

DXDPDT is the matrix containing the time derivatives, — / 

NEQ is the number of columns in ^ ^fp^ matrix. Note NEQ includes 

the derivatives corresponding to the initial 

condition errors. 


6.11 SUBROUTINE DERVX1 

This subroutine evaluates time derivatives of the aircraft state vector 
with respect to the initial condition errors. 


Calling Sequence - CALL DERVX1 (DXDE , DXDED , NX) 


where 


DXDE is the partial derivative, -~— 


DXDED is the time derivative of the above. 


NX is the row dimension of 


m 


matrix. 


d_ /a£\ 
’ dt \3e ) 


Equation: 


dt y 3e J y3e j 
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6.12 SUBROUTINE DERVX2 


This subroutine computes the partial derivatives of the simulated air- 
craft state vector with respect to the control error sources. 

Calling Sequence - CALL DERVX2 (DXH ,DXHD,NX,DUDE) 



6.13 SUBROUTINE DERVY1 

This subroutine evaluates time derivatives of the lagged output measure- 
ment vector with the lag error sources . 

Calling Sequence - CALL DERVY1(DYDE,DYDED,XX,YYL,YYI,UU,NX) 


where 

DYDE 

DYDED 


is the partial derivative, 
is the time derivative of 


9e 

8y L 


9e 


XX is the aircraft state vector, x 
YYL is the lagged output vector, 
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YYI is the indicated output vector, y 
UU is the input control vector, u 

NX is the number of components in output measurement vector, y 


Equation: 


d_(^L 

dt \3e 


)- 


3F 

Jr <Hx + Du - V 


F (- 



6.14 DMTMY1 

This subroutine computes the product of a double precision matrix and a 
single precision matrix and stores the result in a single precision matrix, 

C = A x B. 

Calling Sequence - CALL DMTMY1(A,B,C,L,M,N,N1,N2,N3) 


where 


A is a double precision matrix of dimension N1 by M 

B is a single precision matrix of dimension N2 by N 

C is a single precision matrix of dimension N3 by N 

L is the number of rows of K used to compute the product 

M is the number of columns of A and the number of rows of B used 
to compute the product 

N is the number of columns of B used to compute the product 
N1 is the row dimension of A 

N2 is the row dimension of B 

N3 is the row dimension of C 
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Note that the multiplication is performed only by using the L by M sub- 
matrix of A and M by N sub-matrix of B. Thus, only L by N sub-matrix of 
C is computed. 


6.15 SUBROUTINE ENSMB 

This subroutine is used'to evaluate the ensemble equations for both the 
Ensemble Analysis and the Simulated Data Analysis. 

If the Ensemble Analysis is being performed, this subroutine also computes 
quantities necessary to do the error sensitivity study during the trajectory 
calculation. 

The output from the Ensemble Analysis is then the sensitivity derivative, 

ft)* <“>(!! 

If the Simulated Data Analysis is being performed, then the nominal tra- 

/\ 

jectory data such as the state vector, x, the partial derivatives of the state 

A 

9 x ^ 

vector, , the simulated states, y, and the weighted partial derivatives 

Op 

-i(*A 

of the simulated states, R \8p I , are computed and stored on a scratch file, IU1 
to be used during the Simulated Data Analysis. 


\ 3e / 


and the covariance matrix 



6.16 SUBROUTINE ENSMBL 


This subroutine evaluates ensemble equations by integrating for one 

a 2 . 

sample period and forms I — ^ | matrix. 


( 0 ) 


Calling Sequence - CALL ENSMBL (KKQS, STATES, NSTEP,DELT) 


where 

KKQS is the number of vector equations to be integrated 
STATES is the vector equations 


NSTEP is the number of integration steps until sample time 
DELT is the integration step-size 

6.17 SUBROUTINE ERRSET 

This subroutine sets up a vector of flags ERRNDX corresponding to every 
error source. 

ERRNDX(i) = N =0: means the i^ error source is absent 

4 0: means the i^ error source occupies n C ^ position 
in ERRVEC vector. 

This subroutine also sets up a vector ERRVEC and NAMERR where non-zero 

error value and the names of non-zero error sources are stored. ERRVEC and 

til 

NAMERR are packed arrays. If ERRNDX (i) = N ^ 0, then the i error value is 
stored in ERRVEC (n) and hollerith name is stored in NAMERR (n, 1) and NAMERR (n,2). 

Calling Sequence - CALL ERRSET (BIAS, EP,GM,EIC,FF,NAMPAR,NAMICC) 

where 

BIAS is a vector of bias errors (b's) 

EP is a vector of scaling errors (e’s) 

GM is a vector of misalignment (y's) 

EIC is a vector of initial condition errors (x(0)) 

FF is a vector of lag errors (f's) 

NAMPAR is a vector of hollerith names for parameters 


NAMICC 


is a vector of hollerith names for initial conditions 


6.18 SUBROUTINE ESTNAM 


This subroutine sets up a vector of hollerith names for all parameters 
to be estimated. 


Calling Sequence - CALL ESTNAM 


6.19 SUBROUTINE GET 

This subroutine brings the simulated state vector, measurement vector and 
partial derivatives into memory from the scratch file (IU1) to be used during 
the Simulated Data Analysis. 

Calling Sequence - CALL GET(ITM) 

where 

ITM is the sample data point number. 

6.20 SUBROUTINE GETU 

This subroutine retrieves control vector, u , and the partial derivatives 

with respect to error source,! -g^—1, from the scratch file when control lags 
are present. 

Calling Sequence - CALL GETU (ITM, NS, DUDE) 
where 

ITM is the sample point number 

NS is the number of integration steps to the next sample point 

( 9u l 

DUDE is the partial derivative, 



6.21 SUBROUTINE GETY 


This subroutine retrieves output measurement vector y from the scratch 
file during the Simulated Data Analysis if the control error is present. 

Calling Sequence - CALL GETY(ITM) 

where 

ITM is the sample point number 

6.22 SUBROUTINE INPUTX 

This subroutine controls most of SCIP2 input. This subroutine reads all 
namelist data and constructs F, G, H, and D matrices. 

Calling Sequence - CALL INPUTX 


6.23 SUBROUTINE INTEG 

This subroutine controls the integration of aircraft dynamic equations. 
It integrates one time step, At g , for each call to INTEG. 

Calling Sequence - CALL INTEG (KKQS , STATES , I) 

where 

KKQS is the number of vector equations to be integrated 
STATES is the KKQS vector equations 
I is the sample time point number 

6.24 SUBROUTINE MATADD 


This subroutine performs matrix addition or matrix subtraction. 


Calling Sequence - CALL MATADD (A,B,C,N,M,NA,NB,NC,L) 


where 

A is a NA by M matrix 

B is a NB by M matrix 

C is a NC by M matrix 

N and M specify the submatrix for which addition or subtraction 

is to be performed 

NA, NB, and NC specify the row number of matrices A, B, and C, 
respectively . 

L is a flag to specify addition or subtraction 
=0: addition 
#0: subtraction 

6.25 SUBROUTINE MATMPY 

This subroutine multiplies two matrices, C = A • B. 

Calling Sequence - CALL MATMPY (A, B,C,L,M,N,N1,N2,N3) 
where 

A is a N1 by M matrix 

B is a N2 by N matrix 

C is a N3 by N product matrix 

L is the number of rows in A for which multiplication is performed 

M is the number of columns in A and the number of rows in B for 

which multiplication is performed 
N is the number of columns in B and C for which multiplication 
is performed 

Nl, N2, N3 are the row dimensions of A, B, and C matrices, 
respectively. 


6.26 SUBROUTINE MATPNT 


This subroutine prints out covariance matrices with the parameter names 
top to bottom and across the page. It prints 10 columns at a time. 

Calling Sequence - CALL MATPNT (MATRIX, Nl, N2 , NX, HEAD) 

where 

MATRIX is the covariance matrix 

Nl and N2 specify the dimensions of submatrix to be printed 
NX is the row dimension of MATRX 

HEAD is the header vector. It may contain up to 40 characters 

6.27 SUBROUTINE MEANER 

This subroutine constructs T-matrix for each of mean error sources or 
F m matrix if output lags are present. 

Calling Sequence - CALL MEANER (LERR, IERR) 

where 

LERR is the lag error vector 

IERR is the flag to indicate the end of mean error sources 
=0: continue with mean error sources 
^0: finished mean error sources 


6.28 SUBROUTINE OPDB 

This subroutine controls the print out of input summary for the SCIP2 
program. 

Calling Sequence - CALL OPDB(BIAS,EPS,GM,F,W,EIC,NAMPAR) 
where 

BIAS is the bias error vector 

EPS is the output scale error vector 

GM is the misalignment error vector 
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F 


is the output measurement lag error vector 
W is the output measurement noise vector 
EIC is the initial condition error vector 

NAMPAR is the array containing the hollerith names of the output 
# measurement states 


6.29 SUBROUTINE OTHERR 


This subroutine constructs T-matrix for random error sources. 


Calling Sequence - CALL OTHERR (BERR, EERR.GERR, ICERR, IERTN) 


where 

BERR is the bias error vector 

EERR is the scaling error vector 

GERR is the misalignment vector 

ICERR is the initial condition error vector 

IERTN is the flag to indicate the end of random error sources 
=0: not yet end 

=1: end of random error sources 


6.30 SUBROUTINE PACK 

This subroutine moves vectors into temporary storage area to form a 
packed vector for simultaneous integration. 


Calling Sequence - CALL PACK(V1, IS ,NUM) 


where 

VI is the vector to be moved to temporary storage area 
IS is the starting location -1 in the temporary storage area 
of VI 

NUM is the number of elements in VI to be moved to the temporary 
storage area 
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6.31 SUBROUTINE PACK2 


This subroutine moves two dimensional arrays into temporary storage 
area to form a packed vector. 

Calling Sequence - CALL PACK2(A, J2,N1,N2,NX) 

where 

A is a two dimensional array to be moved to temporary storage area 
J2 is the starting location -1 in the temporary storage area 
of A 

N1 is the number of rows of A to be moved 
N2 is the number of columns of A to be moved 
NX is the row dimension of A 

6.32 SUBROUTINE PERTRB 

This subroutine computes the T-matrix by calling CNSTRT or TCONST for 
a given error source. 

Calling Sequence - CALL PERTRB (ERROR, IB, IERTN) 
where 

ERROR is the standard deviation of the error source 
IB is the error number 

IERTN is the flag to indicate the presence (IERTN=0) or absence 
(IERTN=0) of error number IB 

6.33 SUBROUTINE PRINTF 

This subroutine prints matrix for Ensemble Analysis and (<5p) vectors 

corresponding to random error sources for M0DE=1 of the Simulated Data 
Analys is . 
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Calling Sequence - CALL PRINTF (HEAD, MATRIX) 


where 


HEAD is a header vector of up to 40 characters 

MATRIX is the matrix containing either matrix or (6p) vectors, 


6.34 SUBROUTINE PRINTE 

# This subroutine controls the printout of normalized covariance matrix. 

Calling Sequence - CALL PRINTE (COV, HEAD) 

# where 

COV is the covariance matrix 
HEAD is the header vector 


6.35 SUBROUTINE PRINTS 

This subroutine prints measurement vector and control vector at each 
sample time point. 

Calling Sequence - CALL PRINTS (NN,Y) 


where 

NN is the sample point number 
Y is the measurement vector 


6.36 SUBROUTINE PRINTV 

This subroutine prints the expected values of Sp vector. 

Calling Sequence - CALL PRINTV (HEAD, EVDP) 

where 

HEAD is the header vector 

EVDP is the expected value vector 
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6.37 SUBROUTINE PRTLS1 


This subroutine controls the computation of 


■1^- term for the error 
3e 


sensitivity study of the Ensemble Analysis, 



™ If - If <Hx 


+ Du) + ^ 


9B 

9e 


Calling Sequence - CALL PRTLSl(INDX) 
where 

INDX is the sample point number 


6.38 SUBROUTINE PUT 

This subroutine stores the simulated state vector, the measurement vector 
and their partial derivatives with respect to parameters on a scratch file 
to be used during the Simulated Data Analysis. 

Calling Sequence - CALL PUT(ITM) 

where 

ITM is the sample point number 


6.39 SUBROUTINE RN 

This function computes a normally distributed random number with zero 
mean and unit variance. 

Calling Sequence - RAND=RN(I) 

where 

RAND is the random number returned by the function 
I is a dummy variable (CDC-6600) 


6 . 40 SUBROUTINE RESTOR 


This subroutine either saves or restores the input errors during the 
MODE 2 of the Simulated Data Analysis. 

Calling Sequence - CALL RESTOR(BIAS,EPS,GM,EIC,N) 

where 

BIAS is a bias error vector 
EPS is a scaling error vector 
GM is a misalignment vector 
EIC is an initial condition error vector 
N is a flag to save or restore data 
=0: save 
>0: restore 


6.41 SUBROUTINE RK4 

This subroutine integrates a matrix of differential equations one step 
using a 4th order Runge-Kutta method. 

Calling Sequence - CALL RK4(X,T,DT,K,L,M) 


where 


X is a matrix of the initial values and temporary storage area. 
Upon completion of this subroutine, the integrated matrix is 
stored in place of the initial values. 

T is the integration variable 
DT is the integration step-size 

K is the number of rows of X-matrix to be integrated 
L is 4 times the column dimension of X 

M is the maximum row dimension of X 
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6.42 SUBROUTINE RUNKUT 


This subroutine integrates a vector equation by using a 4th order Runge- 
Kutta scheme. 

Calling Sequence - CALL RUNKUT (N,T,DT,X,DX,KUTTA, SAVD) 
where 

N is the number of elements in the vector 

T is the integration variable 

DT is the integration step-size 

X is the vector function 

DX is the derivative of X 

KUTTA is the counter 

SAVD is the temporary storage area 

6.43 PROGRAM SCIP2 

This is the main program of the SCIP2 program. It controls the input, 
output of input summary, the Ensemble Analysis and the Simulated Data 
Analysis. 

6.44 SUBROUTINE SETNDX 

This subroutine sets up the error index for each error source and the 
error name for the output header. 

Calling Sequence - CALL SETNDX ( A , M, NE , NAME) 

where 

A is the array containing the error source 
M is the dimension of A 

NE is the number of non-zero elements in A 

NAME is the array containing the hollerith names of the error 
source in A 
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6.45 SUBROUTINE SMDATA 


This is the main driver for the Simulated Data Analysis. It calls 
SM0DE1 for M0DE=1 option or SM0DE2 for M0DE=2 option. 

Calling Sequence - CALL SMDATA(EPS , EIC, W, BIAS, GM,FF) 

where 

EPS is the output scale error vector 

EIC is the initial condition error vector 

W is the noise vector 

BIAS is the bias error vector 

GM is the misalignment vector 

FT is the lag error vector 

6.46 SUBROUTINE SM0DE1 

This subroutine performs M0DE1 of the Simulated Data Analysis. 
Calling Sequence - CALL SM0DE1(EPS, EIC, W, BIAS, GM,FF,EPSC,BIASC,FFC) 


where 

EPS is the output scale error vector 

EIC is the initial condition error vector 

W is the noise vector 

BIAS is the bias error vector 

GM is the misalignment vector 

FF is the output lag error vector 

EPSC is the control scale error vector 

BIASC is the control bias error vector 

FFC is the control lag error vector 
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6.47 SUBROUTINE S MODE 2 


This subroutine performs M0DE2 of the Simulated Data Analysis. 


Calling Sequence - CALL SMODE2(EPS,EIC,W,BIAS,GM,FF,EPSC,BIASC,FFC,WC) 


where 

EPS is the output scale error vector 

EIC is the initial condition error vector 

W is the output noise vector 

BIAS is the output bias error vector 

GM is the misalignment vector 

FF is the output lag error vector 

EPSC is the control scale error vector 

BIASC is the control bias error vector 

FFC is the control lag error vector 

WC is the control noise vector 

6.48 SUBROUTINE SYMVRT 

This subroutine inverts a variably dimensional symmetric matrix. The 
calculation is done in double precision. 

Calling Sequence - CALL SYMVRT (A , N, M, DET) 

where 

A is the matrix to be inverted. Upon exit from this subroutine, 
A contains the inverse. 

N is the dimension of submatrix of A to be inverted 

M is the max imum row dimension of A 

DET is the absolute value of the determinant of A 


6.49 SUBROUTINE TCONST 


This subroutine constructs T-matrix, the scaling matrix for the output 
measurements, corresponding to mean errors. 

Calling Sequence - CALL TCONST (INDEX) 

where 

INDEX is the index for the error sources 

6.50 SUBROUTINE TUN IT 

This subroutine constructs T-matrix corresponding to zero scaling error. 
Calling Sequence - CALL TUNIT 


6.51 SUBROUTINE UNPACK 

This subroutine unpacks vectors from temporary storage area. 

Calling Sequence - CALL UNPACK(V1, IS,NUM) 
where 

VI is the vector where values are to be moved from the temporary 
storage area 

IS is the starting location -1 in the temporary storage area where 
VI is stored 

NUM is the number of values to be moved 

6.52 SUBROUTINE UNPAK2 

This subroutine unpacks two dimensional array from the temporary 


storage area. 



Calling Sequence - CALL UNPAK2(A, J2,N1,N2,NX) 


where 


A is the two dimensional array where the values are to be stored 
J2 is the starting location -1 in the temporary storage where the 
values start 

N1 is the row number of the submatrix of A to be filled 
N2 is the column number of the submatrix of A to be filled 
NX is the row dimension of the matrix A 


6.53 SUBROUTINE YPRTLS 


This subroutine computes the simulated state vector, y , 

xv /\ 

3y 3y 

derivatives, .with respect to the parameters, and R 

3 2 J . . 

— -z matrix. 

3p 


its partial 
to form 


Calling Sequence - CALL YPRTLS 


LABELED COMMON DESCRIPTIONS 


VII 


LABELED COMMON DESCRIPTIONS 

The labeled commons used in the program and the variables in the labeled 
commons are described in this section. The BLOCK DATA subprogram includes 
all the labeled commons and defines the length of blocks. The labeled commons 
in any other subroutine may be cut off at the last variable referenced in 
that subroutine. Thus, they may be shorter. 


# 
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Common /BIASES/ B LONG (7) , BLAT(7) 


BLONG 


BLAT 


The bias error for the measured state vector of the longi- 
tudinal analysis; 

BLONG (1) = B 0 

BLONG (2) = B q 

BLONG (3) = B 

a 

BLONG (4) = B u 
BLONG (5) = B 

nx 

BLONG (6) = B nz 
BLONG (7) - B* . 

The bias error for the measured state vector or the lateral 
analysis; 

BLAT (1) - Bg 

BLAT(2) = B 

P 

BLAT(3) = B r 
BLAT (4) = B^ 

BLAT(5) - B ny 
BLAT (6) = B* 

BLAT (7 ) = B* . 


Common /CGNCRT/ ECGVEC(3) 


ECGVEC 


The vector of the center of gravity position error 
ECGVEC (1) = e 

cgx 

ECGVEC (2) = e 

cgy 

cyz ’ 


ECGVEC (3) = e 


Common /CNFLCT/ 
SPLCNF 

FULCNF 

LATCNF 

SPLBCN 

FULBCN 

LATBCN 


SPLCNF (3), FULCNF (4) , LATCNF (4) , SPLBCN(5), FULBCN(7), 
LATBCN (7) 


A conflict vector for the ESTIM=3 short period analysis. 

If an element is set to TRUE, the corresponding initial 
condition error is not estimated but treated as error 
source. 

A conflict vector similar to SPLCNF but for the full longi- 
tudinal analysis. 

A conflict vector similar to SPLCNF but for the lateral 
analysis . 

A conflict vector for the ESTIM=4 short period analysis. 

If an element of SPLBCN is set to TRUE, the corresponding 
bias error is not estimated but treated as error source. 

A conflict vector similar to SPLBCN but for the full longi- 
tudinal analysis. 

A conflict vector similar to SPLBCN but for the lateral 
analysis . 


Common /CNTDAT/ TC(2), BC(2), FC(2), UL(2), UM(2) 


TC 

BC 

FC 

UL 

UM 


A vector containing the diagonal elements of the input 
control scale matrix. 

A vector containing the bias errors for the input control 
vector . 

A vector containing the diagonal element of the F^-matrix 
for the input control lag. 

A vector containing the lagged input control vector. 

A vector containing the sampled input control vector. 
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Common /CONTRL/ NOP, DELE(300) , DELA(300) , DELR(300) 


NOP 

DELE 

DELA 

DELR 


The number of sample points. 

The vector of control input sequence, 

The vector of control input sequence, 6^ 

The vector of control input sequence, 6^ 


DELE, DELA, and DELR are in degrees. 



Common /CPERRV/ WDE, WDA, WDR, BDE, BDA, BDR, EDE, EDA, EDR, FDE, FDA, FDR 


WDE 

WDA 

TOR 

BDE 

BDA 

BDR 

EDE 

EDA 

EDR 

FDE 

FDA 

FDR 


The standard deviation, W^ , of the noise for the input 
control vector, 6^, of the longitudinal equation. 

Similar to WDE but for the control vector, 6 , of the 
lateral equation. 

Similar to WDE but for the control vector, 6 , of the lateral 
equation. 

The bias error for the control vector, 6 . 

’ e 

The bias error for the control vector, 6 . 

’ a 

The bias error for the control vector, 6 . 

* r 

The scaling error for the control vector, 6 e> 

The scaling error for the control vector, 6 . 

cl 

The scaling error for the control vector, 6 . 

An element of the F„-matrix for the control vector, 6 . 

C e 

An element of the F -matrix for the control vector, 6 . 

C a 

An element of the F -matrix for the control vector, 5 . 

C r 


Common /DDDPBK/ R0WDL0(9), R0WDLA(13), C0LDLA(13) 


The vector of row indices for non-zero elements of 


matrix for the longitudinal analysis. 


ROWDLO(i) = n means 




is the only non-zero element. 


The vector similar to ROWDLO but for the lateral analysis 


The vector of column indices for non-zero elements of 


matrix for the lateral analysis. 

ROWDLA(i) = n ) /3D \ 

( means [ -z — J 

COLDLA(i) = m \ dp i/n,m 


( 8 ) 


is the only non-zero element. 
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Common /DFDPBK/ R0WFL0(9), C0LFL0(9) , R0WFLA(13), C0LFLA(13) 


ROWFLO 


COLFLO 


ROWFLA 

COLFLA 


The vector of row indices for non-zero elements of 
matrix for the longitudinal analysis. 

The vector of column indices for non-zero elements of 
matrix for the longitudinal analysis. 

■g 

matrix for the lateral analysis. 

The vector of column indices for non-zero elements of 
matrix for the lateral analysis. 


Common /DGDPBK/ R0WGL0(9), R0WGLA(13), C0LGLA(13) 


ROWGLO 


The vector of row indices for non-zero elements of 
matrix for longitudinal analysis. 

The vector of row indices for non-zero element of | 
matrix for lateral analysis. 


The vector of column indices for non-zero elements of 
matrix for the lateral analysis. 
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Common /DHDPBK/ ROWHLO(9), COLHLO(9), R0WHLA(13), C0LHLA(13) 


ROWHLO 


COLHLO 


ROWHLA 


COLHLA 


A vector of row indices of non-zero elements of 
matrix for the longitudinal analysis. 

A vector of column indices of non-zero elements of 
matrix for the longitudinal analysis. 






($?) 


A vector of row indices of non-zero element of 
for the lateral analysis. 

A vector of column indices of non-zero elements of 
matrix for the lateral analysis. 


matrix 


(g) 
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Common /DMATRX/ DL0NG(7), DLAT(7,2) 



The D-matrix for the longitudinal equation; 



The D-matrix for the lateral equation; 
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Common /ERRMAT/ NERR, THMAT(7,4), PDYDE (7, 45), PDXDE(4,45), PDPDE (24, 45) , 
ERRVEC(45) , ERRNDX (45) 


NERR 

THMAT 

PDYDE 


PDXDE 

PDPDE 


ERRVEC 

ERRNDX 


Number of non-zero error sources . 

The matrix resulting from multiplying the scale matrix 
T and the H-matrix. 

The partial derivatives of the measured states with 
respect to the error sources. Each column of PDYDE 
matrix represents derivatives of state with respect to 
a particular error source. 

The partial derivatives of the aircraft states with 
respect to error sources. 

The sensitivity of parameter estimates with respect to 
error sources. Each column of PDPDE matrix represents 
sensitivity of estimates with respect to a particular 
error source. 

The vector of non-zero errors. All non-zero errors are 
packed into the first NERR location. 

The vector of error source flags. 

ERRNDX (i) = n = 0: means that the i-th error source 

is absent 

^0: means that the i-th error source 

occupies n-th position in ERRVEC 
array. 


> 


Common /EXP MAT/ EVDP(24), EDPDPN(24,24) , PDJDP(24) 


EVDP The vector of expected value of 6p, E(6p) . 

T 

EDPDPN The covariance matrix, E(6p6p ). 

PDJDP The partial derivative of the performance function, 

J, with respect to parameters. 

J *2 <y i - (y i - y ±> 


jhJ 

3P 


ay i\ T -1 

>7 R (y -V 


87 



Common /FMATRX/ FFL0NG(4,4), FFLAT(4, 4) 


FFLONG 


FFLAT 


The F-matrix for the longitudinal equations; 




0 10 
0 Mq Mw 


2 V 

“■ sin 0 ~r~ cos a Zw 

k ok o 


e V 

f cos 0 — i— sin a Xw 

k o k o 


0 

Mu 

Zu 

Xu 


The F-matrix for the lateral equations; 




N, 


sin a 


N 


- cos a 


N 

r 

tan 0 


cos 0 
V o 

0 

0 

0 
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Common /GENERR/ T(7, 7) , TT(7,7), NOT, NDX(7) 


T 

TT 


The optional input scale matrix. 

The scale matrix; for the longitudinal equations, the non 
zero elements are: 

TT(1,1) - 1 + e Q 

TT(2,2) - 1 + e 


TT{3,2) = - 


(e + e ' 
£ gx 


TT(3,3) - 1 + e 
TT (4,4) = 1 + e 


a 


TT(5,5) * 1 + e 
Y 

TT(5,6) —■ 


nx 


TT(5 ,7) » 


e + e 
az cgz 
Kg 


TT(6,5) « 


nz 


TT(6,6) - 1 + e 


TT(6,7) = -I 


nz 

fe + e \ 
ax cg x 




TT(7, 7) = 1 + e* 

for the lateral equations, the non-zero elements are: 
TT(1,1) = 1 + e g 

TT(1,3) - . (!g-1 3s s -) 

TT(2,2) = 1 + e 
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Common /GENERR/ (Continued) 


NOT 


TT(2,3) = - -E. 

Y r 

TT(3, 2) = ~ 

K 

TT(3,3) = 1 + e 

i 

TT (4 , 4 ) = 1 + e 
TT(5,5) = 1 + 




ny 


TT(5, 6) = ^ - X - 

; £ cez \ 

c g / 


TT(5, 7) = - 


TT (6 , 6 ) = 1 + e * 
P 

Y* 

TT (6 , 7) = - ^ 

Y* 

TT(7,6) = 

TT(7,7) = 1 + e* 


A flag to indicate whether or not the error source is present 
in TT matrix. 

=TRUE: error source is present in TT matrix 

=FALSE: no error source present in TT matrix 


NDX 


A vector of indices of the non-zero element of TT matrix. 


Common /GMATRX/ GLONG (4) , GLAT(4,2), XXI, XZI, ZZI 


GLONG 


GLAT 


XXI 

XZI 

ZZI 


The G-matrix for the longitudinal equations 


r on 



e 


L o J 

The G-matrix for the lateral equations; 



L 0 0 


Moment of inertia about x, I 

xx 

Product of inertia about x-z, I 

’ xz 

Moment of inertia about z, I 

’ zz 



Common /HEAD/ 


ENSMBL 


SIMDAT 


LONGNS 


ENSMBL (2) , SIMDAT(3), LONGNS(7), LATNS(7) 


TVo word array containing hollerith data; 

ENSMBL (1) = ' ENSEMB' 

ENSMBL (2) = 'LE** ' 

for print out. 

Three word array containing hollerith data; 

SIMDAT (1) = ’SIMULA' 

SIMDAT (2) = ’TED DA’ 

SIMDAT(3) = ’TA ’ 

Seven word array containing the longitudinai instrument 
names. 

Seven word array containing the lateral instrument names 


LATNS 


Common /HMATRX/ HLONG(7,4), HLAT(7,4) 


HLONG 


HLAT 


The H-matrix for the longitudinal equation; 


0 

0 


0 

1 


M 


0 0 
k cos a -k sin a 


V 


w 


w 


0 — 


M 


w 


u 


0 — — 
8 g 

Z 

g 

M 


u 


The H-matrix for the lateral equation; 


1 

0 

0 

0 

VY 

kg 

L 


3 


6 


L N 6 


0 

1 

0 

0 


0 

0 

1 

0 


0 0 


N 


O' 

0 

0 

1 

0 

0 


N 0 
r 
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Common 


/ HNDMAT / 


H(7 ,4 ) , D(7,4), F(4,4), G(4,2) 


H 

D 

F 

G 


The H-matrix for the case. It may either be HLONG or HLAT 

depending on the input option. 

The D-matrix for the case. 

The F-matrix for the case. 

The G-matrix for the case. 



Common /ICERRS/ EICLON(4) , EICLAT(4) 


EICLON 


EICLAT 


Initial condition errors for longitudinal equation; 
EICLON(l) = 6 Q 

EICLON(2) = 6 

q 

EICLON (3) = 6 w 
EICL0N(4) = 6 . 

Initial condition errors for lateral equation; 
EICLAT (1) = 

EICLAT (2) = 6 p 

EICLAT (3) = 6 r 

EICLAT (4) - 6 . 


f 


i 
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Common /INPUTT/ 


INPUTS (30) 


INPUTS 


Input and output options; INPUTS (1) through INPUTS (22) 
controls input and INPUTS (25) through INPUTS (30) controls 
output. INPUTS (23) is a special option. 


Common /LABELS/ 


FLA 


SPL 


LATAN 


< 


FLA(5), SPL (6), LATAN<3) 


Header array for "Full Longitudinal Analysis" option; 

FLA(l) = 'FULL L' 

FLA (2) = 'ONGITU' 

FLA(3) = 'DINAL ' 

FLA(4) = 'ANALYS' 

FLA(5) = 'IS '. 

Header array for "Short Period Longitudinal Analysis" option; 

SPL(l) = 'SHORT ' 

SPL(2) = 'PERIOD' 

SPL (3) = ' LONG I ' 

SPL (4) = 'TUDINA' 

SPL(5) = 'L ANAL' 

SPL(6) = 'YSIS '. 

Header array for "Lateral Analysis" option; 


« 


LATAN (1) = 'LATERA' 
LATAN (2) = 'L ANAL' 
LATAN (3) = 'YSIS ' 



Common /LAGERR/ FLONG(7) , FLAT (7) 


FLONG 


FLAT 


The diagonal elements of F^-matrix when output measurement 

lags are present in the longitudinal equations. FLONG 

contains f Q , f , f , f , f , f , and f*. respectively. 

0 q a* u’ nx* nz* q 

The diagonal elements of F^-matrix when output measurement 

lags are present in the lateral equations. FLAT contains 

f a , f , f , f,, f , f*. and f*, respectively. 

0’ p* r* <p ny* p*. r* v J 


Common /LIMITS/ JEQS, KEQS, LEQS, MEQS, NEQS, NUM, MX1, MX2, NSTNDX(7), 

MONTE, NUMERR(II) 


KEQS 


LEQS 


Temporary variable, NEQS + 1. 

9 X 

Number of vector equations, , to be integrated in the 
ensemble equations 

Temporary variable, NEQS + MX2. 


MEQS 


Total number of parameters, including initial conditions 
and bias errors to be estimated by the program. 

Number of parameters to be estimated when ESTIM = 1. 

Number of output measurement instruments. 

Number of elements in the state vector, x. 


NSTNDX 


MONTE 


NUMERR 


A vector of flag to indicate which instruments are present 

NSTNDX(i) = 0: indicates that i^ instrument is absent 

NSTNDX (i) = j > 0: indicates that i^ instrument is the j 

instrument among those present. Note: i j . 

Flag to indicate whether dynamic equations are being integrated 
by the ensemble analysis program (if false) or by the simulated 
data analysis program (if true). 

A vector containing the number of error sources; 

NUMERR (1) = NB Number of non-zero bias error sources 
NUMERR(2) = NE Number of output scale error sources 

NUMERR (3) = NG Number of misalignment error sources 

NUMERR (4) = NF Number of output lags 
NUMERR(5) = NCGA Number of error sources in center of 
gravity position 
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i 


Common /LIMITS/ (Continued) 


NUMERR (6) = NCGB 

NUMERR(7) = NIC 
NUMERR(8) = NBC 

NUMERR(9) = NEC 

NUMERR (10)= NFC 
NUMERR(11)= NCR 


Sum of accelerometer position error sources 
and the vane position error source* 

Number of initial condition error sources. 
Number of bias error sources for input 
control. 

Number of scale error source for input 
control. 

Number of input control lag error sources. 
Number of random error sources for input 
control. 
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Common /LOCCG/ 


EAX, EAY, EAZ, EVX 


EAX \ 

EAY J. Accelerometer location error for x, y and z component 

EAZ ) 


EVX 


Vane location error 


Common /LOGICL/ 


FULTST, SPLTST, LATEST, HERAT 


FULTST 


SPLTST 


LATEST 


ITERAT 


Logical variable, if true, to indicate that the full 
longitudinal equations be executed. 

Logical variable, if true, to indicate that the short 
period longitudinal equations be executed. 

Logical variable, if true, to indicate that the lateral 
equations be executed. 

Logical variable, if true, to indicate the SIMULATED DATA 
phase of the INPUTS (23) option. 


Common /MEANIT/ IA2, IB2, IC2 


IA2 
IB 2 
IC2 


Pointer to run down the input control error sources. 
Pointer used to pick out a particular input control error. 
Not used. 


« 


« 
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Common /MISALN/ GMLONG(7), GMLAT (7) 

GMLONG Misalignment error sources for the longitudinal equations; 

GMLONG(l) = y 0 

GMLONG (2) = Y q 

GMLONG (3) = Y a 

GMLONG (4) = Y u 

GMLONG (5) = Y 

nx 

GMLONG (6) = Y 

nz 

GMLONG (7) « Y* • 

q 

GMLAT Misalignment error sources for the lateral equations; 

GMLAT (1) = Yg 

GMLAT(2) = Y 

P 

GMLAT (3) = Y r 
GMLAT (4) = y^ 

GMLAT (5) - Y ny 
GMLAT (6) = Yp 
GMLAT(7) = Y*- 
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Common /MISCEL/ DT, NMC, GRAV, UNITS, RADDEG, NDXP, PERT, TIME, DTS , TAU 


DT 

NMC 

GRAV 

UNITS 

RADDEG 

NDXP 


PERT 

TIME 

DTS 


Integration step-size. 

Number of Monte Carlo runs to be executed in the Simulated 
Data Study. 

Acceleration due to the gravity. 

Hollerith data for the unit of the gravity. 

Conversion factor from degrees to radians. 

Flag used to perturb parameters under INPUTS(23) option 
for the lateral equations. 

=0: perturb all parameters 

1 £ n _< 13: n fc ^ parameter by the amount specified by PERT 

Amount in percent the parameters are to be perturbed. 

Time of the trajectory. 

Sample rate . 

Time to the first sample point of the input control variable. 


TAU 


Common /MISCL1/ NUM1, NUM2, MODM2, II, SWICH, IU1, IU2 


NUM1 

NUM2 

MODM2 

II 

SWICH 

IU1 


Number of components in the control variable, 
=1: for the longitudinal equations 
=2: for the lateral equations 

Not used. 

Variable MODE-2. 


Running counter of the Monte Carlo runs. 


Flag, if non-zero, to indicate that the initial condition 
error is being considered. 


Scratch unit to save x, -tt - 

8p 


r - 1 9Z 
3p 


IU2 


Scratch unit to save 6p for each of M0DE=2 Monte Carlo 


Common /NAMVAR/ 


NAMVAR(7,2), NAMIC(4,2), NAMVAN(4), NAMCG(3) , NAMERR(45, 2) , 
NAMEST(24) 


NAMVAR 

NAMIC 

NAMVAN 

NAMCG 

NAMERR 

NAMES T 


An array containing the hollerith names of measurement 
states for the longitudinal and the lateral equations. 


An array containing the hollerith names of dynamic states 
for the longitudinal and the lateral equations. 


An array containing the hollerith names of the accelerometer 
location and vane location error sources. 


An array containing the hollerith names of the center of 
gravity location error. 

An array containing the hollerith names of the error 
sources. 

An array containing the hollerith names of the parameters 
being estimated. 
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Common /NOISE/ 


WLONG 


WLAT 


WLONG (7) , WLAT (7) 


An array containing the standard deviations of the white 
measurement noise for the longitudinal states. 

An array containing the standard deviations of the white 
measurement noise for the lateral states. 



Common /ESTATE/ XN(4), PDXDP(4,17), YE(7), PDYDP(7,24) 


XN 

PDXDP 


YN 

PDYDP 


A 

The simulated aircraft state vector, x. 


The partial derivatives of the simulated aircraft state 
with respect to parameters and the initial conditions, 


3x 

9p 


A 

The simulated measurement state, y. 


The partial derivatives of the 
with respect to parameters and 


simulated measured state 
the initial conditions, 


li 

9p * 
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Common /OPTNS/ ENORSD, EQNS, ESTIM, MODE, INSLON(7), INSLAT(7) 


ENORSD 

EQNS 

ESTIM 


MODE 

INSLON 

INSLAT 


Option word containing either Hollerith data ENSMBL or 
SIMDAT depending on whether the Ensemble Analysis or the 
Simulated Data Analysis is to be exercised. 

Option word containing either hollerith data FULL, SPL, or 
LAT depending on whether the Full Longitudinal equations, 
the Short Period Longitudinal equations, or the Lateral 
equations are to be used. 

Option to specify the parameters to be estimated. 

=1: parameters only 

=2: parameters and the initial condition errors. 

=3: parameter, bias errors, and those initial conditions not 
correlated with bias errors 

=4: parameters, initial condition errors, and those bias 
errors not correlated with initial condition errors. 

Option specifying the mode of the Simulated Data Analysis . 

An input vector of the instrument names of the longitudinal 
measured state . 

An input vector of the instrument names of the lateral 
measured state. 


Common /OTHRIT/ IAl, IB1, IC1 


IA1 

IB1 


Variable used to run down the output measurement error 
sources . 

Variable used to pick out an error source. 


IC1 


Variable used to run down the non-zero error sources. 


Common /OUTDAT/ FM(7), BVEC(7), WVEC(7), YL(7) 


FM 

BVEC 

WVEC 

YL 


A vector containing the diagonal element of the F 
for the output measurement lag. 

A vector containing the random bias errors. 

A vector containing the random noise errors. 

A vector containing the lagged output state. 


. -matrix 
M 


Common /REFTRJ/ THETAO, ALPHAO, VEL 


THETAO 

ALPHAO 


Initial pitch altitude, 0 
Initial angle-of-attack, 


VEL 


Total velocity, V. 


Common /SCALN G/ EPL0NG(7), EPLAT ( 7 ) 


EPLONG 


EPLAT 


€ 


Scale factors for the longitudinal measurement; 
EPLONG (1) = e 3 
EPLONG (2) = e^ 

EPLONG (3) = e a 
EPLONG (4) = e u 
EPLONG (5) = e 

nx 

EPLONG (6) = e 

nz 

EPLONG (7) = e* . 

q 

Scale factors for the lateral measurement; 

EPLAT (1) = 

EPLAT (2) = e 

P 

EPLAT (3) = e 

r 

EPLAT (4) = e , 

<P 

EPLAT (5) = e 

ny 

EPLAT (6) = e- 
EPLAT (7) = e* . 
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Common /STATES/ X(4), U(2) , R(7), D2JDP2(24,24) 


X 

U 

R 


The aircraft state vector, x. 

The input control, u. 

The diagonal element of the weighting matrix R ^ , 


R 1 ^ 



The inverse of the second partial derivatives of J with 


respect to the parameters, 



D2JDP2 


APPENDIX A 

SUMMARY OF EQUATIONS CODED IN ANALYSIS PROGRAM 


A. 1 Explanation of Time Increments Used by Program 

The three time steps which are used in the program are: 

At - This is the time step of the Runge-Kutta integration 

Li 

package. It is set small enough so the effect of the 
highest frequency dynamics is correctly simulated. 

This term is usually governed by the control or output 
measurement lag with the smallest time constant. 

At g - This is the sample time increment used by the identifi- 
cation process. When measurement lags are ignored, it 
is also used as the integration step. 

T - This is the time delay that governs when sampling of the 
control input is made after time zero. Normally, 

0 < t £ At . Use of T enables sampling the control input 
at different time points than the output measurements. 

A. 2 Simulated Data Analysis Subprogram - Mode 1 


A. 2.1 Simulated Airplane Equations 

The equations governing the simulated airplane are: 

x = Fx + Gu. ; x(0) = 

. input 0 


(A. 1) 


y = Hx + Du, 

T input 

(true output) 


(A. 2) 
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(indicated output - no lags) 


y L ' 'Vl + Vl ! V 0) ' y i (0) 

(effect of lags) 


(A. 4) 


y i ‘ y Li + *i 

(recorded output. If lags are present, this 

samples y every At seconds starting at 

time At .) s 

s 


A. 2. 2 Control Input Equations 

The equations governing the recorded control input are: 


u 


u input 

(true control which is input to program) 


u = T u + B 
I c c 


(indicated control with no lags) 


(A. 6) 


(A. 7) 


«L = -F c u l + F c u i; u L (0) = (0) 

(effect of lags) 


(A. 8) 


u 


mi 


= u. 


Li 


+ w 


cx 


(A. 9) 


(recorded input. If lags are present, this samples 
u T every At_ seconds starting at time t.} 

L S 


* Note The control measurement noise is only used in the Monte Carlo 
Mode (Mode 2) of the Simulated Data Analysis. It is set 
to zero otherwise (Mode 1) . 



A. 2. 3 Ensemble Analysis Subprogram Equations 


The random errors are first analyzed. To do this, the follow- 
ing equations are first integrated and evaluated in the ensemble 

/ 

analysis subprogram: 


x = Fx + Gu ; x(0) = x Q 

(At integration step) 
s 



(for each initial condition estimated) 





A 


X 


+ 



u 


(from Eq. A. 2) 


3£ „ H 3£ 

3 ?ic 3p ic 

(for each initial condition estimated) 


3y 

3p b 


l 


(for each bias term estimated) 


(A. 10) 


(A. 11) 


(A. 12) 


(A. 13) 


(A. 14) 


(A. 15) 


n r m - 

V #i R -1 #1 

£-> 9p 3p 


(A. 17) 


Invert 


3 J A 


A. 2. 4 Simulated Data Analysis Subprogram Equations for Random Errors 

Next, each of the random output measurement errors is 
sequenced. The object is to evaluate: 


‘2 <y i ' y i )T R_1 


(A. 19) 


(A. 20) 


First, the effect of each unestimated output measurement bias 
is computed. This affects only the residual of Eq. (A. 19) as: 


Cy i - ty ■ y u - y i * 


(A. 21) 


Then, each of the output transformation errors of T is sequenced. 
Again, the effect on Eq. (A. 19) is: . 


( y i - 9 i ) = T(Hx ± + Du.) - (Hx i + Du.) = (T - I) (Hx^ + Du ± ) (A. 22) 
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[ 


Then, each of the unestimated initial condition errors is sequenced. 

This requires the integration of Eq. (A.l) with x(0) = x^. Then, from Eq.(A.2), 


(y ± - y ± ) = (y T± - y ± ) = - j^) 


(A. 23) 


Next, the control measurement errors due to T and B are sequenced. 

, c c 

To do this, the previous results of integrating Eq. (A. 10) are first 
stored. These state values represent the true airplane equations 
in this case. So 


X i " X i 


(A. 24) 


Then, Eqs. (A. 10), (A. 11), and (A. 12) are reintegrated using the At in- 
tegration step and modified u .. Equations (A. 12) - (A. 18) are then 

mi 

reevaluated. For bias errors 


u . = u. + B 
mi i c 


(A. 25) 


For scale factor errors 

u = T u. 
mi c i 


(A. 26) 


For the random errors, the total parameter covariance is 


E Hplp') - E nolse + WPjSPj 1 ) 


(A. 27) 


where r is the number of random error sources. E is the co- 
noise 

variance of the parameter estimate errors due to the output noise w^. 
^noise e * t * ier set to Eq • (A. 18) or computed using the Monte Carlo 
option with 


(y i - ? i> - "i 


(A. 28) 


at each time point. 
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A. 2. 5 Mean Errors 


The mean errors are also sequenced. These include some 

elements in the T matrix (Eq. (A. 3) and all elements in F and F 

m c 

(Eqs. (A. 4) and (A. 8)). If there are no lag errors input, the. evalua- 
tion is just like Eq. (A. 22). 


If there are control measurement lag errors, Eq. (A. 8) needs 

to be evaluated with u = u, over each interval At . The integra- 

X s 

tion step is At^. The stored u^ is sampled every At g but begin- 
ning T seconds after time zero. 


If there are output lags, Eq. (A. 4) needs evaluation. This 
requires reintegration of Eq. (A.l), with x(0) = x^, and Eq. (A. 4) of 
the form 


y = -F y T + F (Hx + Du) 
J L m L m 


(A. 29) 


The step size is At . The output of Eq. (A. 29) is sampled every 

At seconds with 
s 



(A. 30) 


(A. 10) - (A. 18) need to be reevaluated in case of control 


Again, Eqs 
lags. 

however. The effect on Eq. (A. 19) is 


Equations (A. 10) - (A. 12) use an integration step size of At , 

s 


v y i 


y Li ' <H *i + Du »i > 


(A. 31) 


A. 3 Simulated Data Analysis Subprogram - Mode 2 


The random errors contained in B, part of T, B„ T , and x are 
generated at the beginning of each run using the input standard devia- 
tions and the random number generator. They are held constant during 
each single Monte Carlo run. The random noise w^ and w c ^ are regener- 
ated each sample point during each run. Each of the mean errors in T 
plus elements of F^ and F^ are set equal to the values input and are 
not changed during any of the runs. 
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For output measurement errors in T, B and only, the 
residual (y^ - y^) in Eq. (A. 19) is computed by 

y - y ± = 1(11^ + Du i ) + B + w ± - (Hx^ + Du^ (A. 32) 

For random initial conditions added, Eq. (A.l) must be integrated 
each time and Eq. (A. 32) gets changed to 

y i - 9 - T(Hx t + Du i ) + B + w ± - (Hx ± + Du^ (A. 33) 

For non- lag control measurement errors, Eqs. (A. 7) and (A. 9) 
get combined so that at each sample point 


mi 


T u + B 
c c 


+ w 


ci 


(A. 34) 


where w . is randomly generated each time point. Because of this 
cl 

change, Eqs. (A. 10) - (A. 12) require integration each pass through, 
and Eqs. (A. 13) - (A. 18) require re-evaluation each pass through. With 
these changes, Eq. (A. 33) becomes 


y i - ?i - T(Hx i + Du^ + B + WjL - (Hx i + Du mi ) 


(A. 35) 


For lags in the control input measurements and output measure- 
ments, two integration step sizes have to be used. First, Eqs. (A.l) - 
(A. 4) are evaluated using the At T step. The resulting y T is sampled 

L Ju 

# every At seconds starting at time At . Then, Eq. (A. 5) is evaluated 

s s 

and the resulting y stored. Next, Eqs. (A. 7) - (A. 8) are evaluated using 

Lt 

the At. step. The resulting u is sampled every At seconds starting 

L L S 

at time T. Then, Eq. (A. 9) is evaluated and the resulting u . stored. 

ml 

# Next, Eqs. (A. 10)- (A. 12) are evaluated using the At step and the stored 

s 

u .. Then, Eqs. (A. 13) - (A. 18) are evaluated. Then, the residual for 
ml 

Eq. (A. 19) is computed using 

y i ‘ h = y i " + Du mi ) (A. 36) 
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In Mode 2, w . is used lust like w . . 

ci i 


In Mode 2 , several Monte Carlo runs are made. In 

each run, the values of the random elements in T, B, and x are ran- 

domly generated at the beginning of the run. For comparison, the 

standard deviations used to generate these values are the same as 

those values used one at a time in Mode 1. These random values 

remain constant throughout a run, but change from run to run. The 

white noise vector w^ is randomly generated, and it changes every 

sample point throughout the run. Those mean error sources in T, 

F , and F which aren't random are held constant throughout all the 
me 

runs. The error Ap^ in the parameter vector obtained from each run 
is saved. For m Monte Carlo runs, the mean error in p is 


m 


Ap A E {Ap} = ^ J 

j=l 


Ap, 


(A. 37) 


The standard deviation about this mean is 

m 


{6p6p T } = i ^ (Ap - Ap) (Ap - Ap) 


— T 

A \ A 


j=l 


(A. 38) 


For small values of the error sources, the results of Eq. (A. 38) should 
match those of Eq. (A. 27). Also, Eq. (A. 37) should match the value of 
the mean obtained in Mode 1. 
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A . 4 Ensemble Error Analysis Subprogram 


A. 4.1 Output Errors 


The basic equation which is used to compute changes in the 


parameters is 


2 , 1-1 


" 3 J " 
a 2 

L3p J 


r v -i 3 ?i ) i v -i i 

+ ^ Ur * ar I f 1 <vV 

L i=l ) J i=l 


The y are the sampled measurements taken from the aircraft 
(Eqs. (A.l) - (A. 5)) and the ^ are the simulated values obtained from 

the measured control input (Eqs. (A. 6) - (A. 9) and Eq. (A. 10)). The 
term 3y^ comes from Eqs. (A. 11) - (A. 16). 

3p~ 

It is assumed that if lag errors exist, lag terms are present 
in all control input and system output measurement equations . If the 
measurements of y^ have lags, the dynamic equations can be written as 

x = Fx + Gu 


y. = — F y. + F (T [H x + D u] + B ) 
L mi- m 


Otherwise, the output equations are of the form 


y = T(Hx + Du) + B 


(A. 40) 
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The sensitivity of y_^ to each error e is computed. For 


biases , 


£ * 1 


For errors due to the matrix T, the sensitivity can be found as 


£ ■ £ c* + Do) 


For initial condition errors, 


ix 

= h|S 


3e 


3e 


the 

term 

3x 

3e 

comes 

d 

Ax' 

\_ 


dt 

l 3e 


F 3e * 


comes from integrating 
3x 


3e 


(0) = I 


When output lags are present, the sensitivity differential 
equation comes from Eq. (A. 39), and is 



The results of Eq s. (A. 41) through (A. 45) are combined into a general 



3e 


(A. 41) 


(A. 42) 


(A. 43) 


(A. 44) 


(A. 45) 


for each error e which affects the output measurements. 


.4.2 Control Input Errors 


With respect to Eq. (A. 20), the error effect or the term 
is computed. The differential equations of the simulated aircraft are 


’ - F c“l + VV + B c ) ■ U L (0) " U I C0) 


c L c 
x - Fx + Gu Li 


y = Hx. + Dll 
J i 1 Li 


For control measurement lags the following equations, that come 
from Eqs . (A. 46 and (A. 47), are first integrated. 


d_ 

dt 


(§) 

'An I 


3u 

F 3x + G • 
F 3e 3e » 


|| (°) = ° 


9u L 9F c . 3u 

— , ,, . = -F -5 + -r— ' (-U + u) ; — h (0) 

dt \ 3e } c 3e de L ’ a e 


3u^ 
9e . 


(A. 46) 
(A. 47) 
(A. 48) 


( 0 ) 


(A. 49) 
(A. 50) 


Then, from Eq. (A. 48) 




9u Li 

i _ V, i a. n — 

3e H 3e D 3e 


(A. 51) 


For control scale factor errors, the program first integrates 

9T 


d • c 

dt 9e 3e 9e u ’ de 


c 

u ; (0) - 0 


(A. 52) 


Then, the sensitivity is computed as 


9y. 3X 

i. = H * + D 2_c u. 

9e 9e 3e 


(A. 53) 


For control biases, the program integrates 


3_ 

at 



F 
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Then, the sensitivity becomes 


ii 

ae 
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3e 


(A. 54) 


(A. 55) 


A. 4. 3 Effect on Parameter Estimates 


The sensitivity of parameter estimates due to output measure- 
ment is of the form 


3 


(Sp) 


n 



(A. 56) 


3 y, 


where 7 


3e 


comes from evaluating Eqs • (A. 41) through (A. 45). 


The sensitivity of parameter estimates due t:o control input 
errors is of the. form 




3V 1 

La ? 2 


-l n 
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v 

3p 
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-1 


r. A 

th 

3e 


(A. 57) 


where 


Hi. 

9e 


comes from evaluating Eqs. (A. 51), (A, 53), and (A. 55). 


For random errors, the. total parameter covariance is 


E ({p 6p T > - E nolse + -i <{p)E {e e 1 ) A (Sp) 1 

For mean errors, the mean parameter error is 
E { 6p } = E { e } 


(A. 58) 


(A. 59) 
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